On the viability of local criteria for chaos 
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Abstract 



We consider here a recently proposed geometrical criterion for local instability based 
on the geodesic deviation equation. Although such a criterion can be useful in some 
cases, we show here that, in general, it is neither necessary nor sufficient for the 
occurrence of chaos. To this purpose, we introduce a class of chaotic two-dimensional 
systems with Gaussian curvature everywhere positive and, hence, locally stable. 
We show explicitly that chaotic behavior arises from some trajectories that reach 
certain non convex parts of the boundary of the effective Riemannian manifold. Our 
result questions, once more, the viability of local, curvature-based criteria to predict 
■^J - ! chaotic behavior. 
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1 Introduction 



In recent years, many efforts have been devoted to formulate a local geomet- 
rical criterion which could indicate chaotic behavior in a mechanical system. 
The main idea behind such attempts is to provide the Lagrangian or Hamil- 
tonian formulation of the equations of motion with an effective Riemaniann 
structure, and then search for some criterion based on its curvature properties 
which could predict local instabilities that might lead to chaos [1,2]. This line 
of research has been mainly inspired by the celebrated result due to Anosov[3] 
stating that the geodesic flow is chaotic in compact manifolds for which all 
sectional curvatures are everywhere negative. 

In a recent work, Zscz§sny and Dobrowolski[4] rederive in a straightforward 
way the relevant results necessary to formulate such criteria. We adopt here 
their conventions, which we will briefly recall now. Let us consider a classical 
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mechanical system with N degrees of freedom described by the Lagrangian 

£(w) = ^(?)9Y'-n9), (1) 



where i,j = 1,2, ... ,N, the dot stands to differentiation with respect to the 
time t, and gij is the Riemannian metric on the configuration space A4. All 
the quantities here are assumed to be sufficiently smooth. The Euler-Lagrange 
equations of (f) can be written as 

7 • -g'O.Viq). (2) 

where r*- fe is the Levi-Civita connection for the metric gij. 
The Hamiltonian of the system described by (1) 

H(q,p) = ±g ii (q) Pi p j + V(q), (3) 

with pi = Qijqi , is obviously a constant of motion, namely the total energy 
For a fixed energy E, the trajectories in the phase-space are confined to the 
hypersurface E = \g % ^PiPj + V(q). The admissible region for the trajectories 
in the configuration space is, therefore, given by 

V E = {q E M : V(q) < E}. (4) 



In general, the region V E can be bounded or not, connected or not. The 
boundary of the admissible region for the trajectories is given by 

dV E = {q E M : V(q) = E}. (5) 



If the potential has no critical points on the boundary (VV ^ 0), then dV E is 
a (^ — ^-dimensional submanifold of M.. We can easily see that if a trajectory 
reaches the boundary &D E at a point qo, its velocity at this point vanishes and 
the trajectory approach or depart from qo perpendicularly to the boundary 
dV E . In particular, there is no allowed trajectory along the boundary. 

The crucial observation here is that the equations of motion (2) are, in the 
interior of T> E , fully equivalent to the geodesic equation in the "effective" 
Riemannian geometry on M. defined by the Jacobi metric [4] 

g ij (q) = 2(E-V(q))g ij (q). (6) 
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The geodesic equation for the effective geometry defined from (6) is given by 

vti-^ + p^^-o m 



where w = dq l /ds is the tangent vector along the geodesic, V and f*- fc are, 
respectively, the covariant derivative and the Levi-Civita connection for the 
Jacobi metric §ij, and s is a parameter along the geodesic obeying 

% = 2(E-V(q)). (8) 



It can be shown that s(t) is a monotonically increasing function of time in the 
interior of Pe[4]. However, from (8) we can see that s is not a good evolution 
parameter on the boundary, since ds/dt = on 8T>e- As we will see, it is 
namely this failure, on the boundary &De, of the equivalence between the 
equation of motion (2) and the geodesic equation in the effective Riemannian 
manifold (7) that precludes the possibility of obtaining an efficacious local, 
curvature-based criterion to predict chaotic behavior of the system governed 
by (2). 

In a Riemannian geometry, the local tendency of geodesies to converge or to 
diverge is governed by the geodesic deviation equation: 

V u V u n = R(u,n)u, (9) 



where R(u, n) is the Riemannian curvature of the Jacobi metric (6). The vector 
field n, called geodesic deviation, is orthogonal to u and measures the distance 
between nearby geodesies. For two-dimensional systems, one can choose Fermi 
frames (E 1 ,E 2 ) such that u — E 1 and n = xE 2 , and the geodesic deviation 
equation takes the simple form [4] 

d 2 x 

2| = -Kx, (10) 



where K is the Gaussian curvature of the two-dimensional Jacobi metric (6), 
which is given by 

1 VV 1 (W) 2 , x 



where K is the Gaussian curvature of the original two-dimensional configura- 
tion space M. The Laplacian V 2 V and the quadratic term (VV") 2 are defined 
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with respect to the configuration space metric gij. Note that, with the hy- 
pothesis that the potential is smooth (V 2 V is finite) and has no critical points 
(VV 7^ 0) on &D E , the effective Gaussian curvature K diverges on the bound- 
ary. Moreover, since E — V > inside V E , K always assume arbitrary large 
positive values near the boundary. It is clear from (10) that when the Gaus- 
sian curvature is positive the geodesies converge locally. On the other hand, 
if it is negative, one has the local divergence of nearby geodesies. Zscz§sny 
and Dobrowolski illustrated the case of positive Gaussian curvature K with 
three known integrable 2-dimensional systems: quadratic potentials, the Toda 
lattice and the bounded Kepler problem. 

If the admissible region for the trajectories on the configuration space has 
no boundary, the equivalence between (2) and (7) is complete. If, besides, all 
sectional curvatures are everywhere negative, close trajectories governed by (2) 
tend to diverge at every point. However, even in such case, it does not mean 
that the system is chaotic, as Zscz§sny and Dobrowolski have shown with their 
unbounded Kepler problem example. The key requirement of compactness 
is lacking here to reproduce the result of Anosov. It is quite easy, indeed, 
to generate integrable models with effective Gaussian curvature everywhere 
negative on a non-compact V E by exploring harmonic potentials in higher 
dimensions [4]. The local divergence of nearby trajectories is not, therefore, a 
sufficient condition to the occurrence of chaos. 

In this work, we question again the viability of such local criteria by show- 
ing, by means of an explicit example, that they are neither necessary to the 
occurrence of chaos. In the next section, we will introduce a two-dimensional 
system with the Gaussian curvature K everywhere positive in the interior of 
T>e, and, hence, locally stable. As we will see, chaotic behavior is indeed ruled 
out for large sets of trajectories reaching the boundary &D E . However, we 
identify a set of trajectories that reach a non convex part of the boundary and 
give rise to chaotic motion. The situation here is similar to case of classical 
billiards[5]. Segments of trajectories between two successive bounces on the 
boundary &D E are integrable, but the smooth match of an infinite number of 
some of these segments give rises to a chaotic trajectory. As, in general, we 
cannot decide a priori if a given trajectory does or does not reach the non 
convex part of boundary for some t G (— oo, oo), the positivity of K inside V E 
is useless to predict if such trajectory is regular or not. 

We notice that the potential problems with the interpretation of the true 
motion governed by (2) as a geodesic equation in an effective Riemannian 
manifold with boundary were already pointed out in [6,7,8]. Some suggestions 
to avoid them based on distinct differential spaces instead of differential man- 
ifolds were proposed in [9]. We also notice that the power of local criteria 
to predict chaos in General Relativity has also been questioned. Vieira and 
Letelier[10], for instance, showed that the local criterion based on the eigen- 
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values of the Weyl curvature tensor presented by Sota, Suzuki and Maedafll] 
were neither necessary nor sufficient for the prediction of chaos in vacuum 
static axisymmetric spacetimes. 

Despite the fact that such curvature-based criteria are neither necessary nor 
sufficient to the occurrence of chaos, the analysis of the geodesic deviation 
equation for the Jacobi metric has proved to be valuable in some cases. For 
instance, in [12], some chaotic behavior, confirmed by the numerical evaluation 
of Lyapunov exponents, is identified in a Newtonian gravitational problem 
with negative effective Gaussian curvature. In [13], chaotic mixing behavior 
in thermal-equilibrium beams, confirmed also by some numerical calculations, 
is again identified in regions with negative Gaussian curvature. In the other 
hand, regions of regular motion in the anisotropic Manev problem[14] and 
the Yang-Mills-Higgs[15] system coincide with positive Gaussian curvature 
regions. We emphasize, however, that since the criteria based on the geodesic 
deviation equation are neither necessary nor sufficient for the occurrence of 
chaos, they must be used merely as preliminary indicators of integrability or 
chaos in dynamical systems. 



2 The model 



Our model is quite simple. It consists in a point-like particle moving under the 
potential corresponding the superposition of a quadratic attractive potential 
and some two-dimensional harmonic repulsive monopoles. The configuration 
space M. is locally just R 2 with the Euclidean metric. Let us consider explicitly 
the simplest case with only one repulsive monopole with charge q, placed at a 
distance a from the center of the quadratic potential. The particle motion is 
governed by the Hamiltonian H = -^v\ + + V(x, v)-, with 

V{x ,y) ^(^ + 4-| ln p±J<!) (12) 



The parameter r has no dynamical consequences, it is chosen in order to 
guarantee that the minimum of V(x, y) be zero. Clearly, in the limit a^Owe 
have a central-force integrable model. The equations of motion are very simple. 
By introducing some dimensionless quantities (x — > ax, y — > ay, t — > ^Jm/kt), 
one has 



a 



p y = y, Pv = -y[i- 2 , 2 

V x 2 + y l 



(13) 



where a = q/ka 2 . 






(b) - 







Fig. 1. The admissible region S i? 2 for the trajectories: (a) the simply-connected 
case E < Eq and (b) the multiply-connected one E > Eq. The boundaries shown in 
(a) correspond to the cases V/Eq = 1/5, 1/2 and 4/5, while in (b) V/Eq = 1.02. For 
all cases, a = 1/2. The non convex parts of the boundaries, including the interior 
part for the multiply-connected case, are the responsible for the bouncing of chaotic 
trajectories. 



We have two kinds of admissible regions for the trajectories. For E < E c , 
where 

£c = f + aln|±j, (14) 



with (3 = \/l + 4a, T>e is simply-connected, in contrast with the case for 
E > Ec (See Fig. 1). In both cases, V(x,y) is smooth in T>e and has no 
critical points on &De- 

In the interior q{T>e, the Gaussian curvature K for this system, 

K i , ("( 1 -^)- 1 ) 2 + ^ 2 ( 1 -^f) 2 n ^ 

2(E-V(x,y)y + 4(E-V(x,y)r ' 1 ° j 



is everywhere positive. Therefore, the nearby trajectories of our model ap- 
proach each other at all points of the interior of V E - As expected, the Gaussian 
curvature diverges on the boundary. 
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3 The dynamics 



In order to study the phase space of our system, we solve numerically the 
system governed by (13) and construct Poincare's sections by using Henon's 
trick[16]. We could do it very accurately, with a cumulative error, measured 
by the constant Hamiltonian 7i, inferior to 10~ 12 . Our results are shown in 
the next figures. 




Fig. 2. Typical Poincare's section (x,p x ) across the plane y = for the sim- 
ply-connected T>e case. For this section, one has a = 1/2 and E/Ec = 4/5. 

In Fig. (2), we have a typical Poincare's section across the plane y = for 
the the simply-connected V E case. The regular solutions correspond to tra- 
jectories that do not probe the non convex part of the boundary, see Fig. (3). 
Despite our exhaustive simulations, we could not find a single trajectory that 
never touches the boundary OT>e- We were not able to prove analytically the 
absence of such trajectories, but we notice that the argument used in [14] to 
prove the existence of them does not hold here. The behavior of the chaotic 
trajectories bouncing on the non convex part of the boundary resemble classi- 
cal billiards [17]. Two close trajectories reaching a convex part of the boundary 
are bounced in a "focusing" way, while the non convex part defocuses nearby 
bouncing trajectories, causing the divergence of the associate geodesies, in 
spite of the positive Gaussian curvature K. 

Similar results hold for the multiply-connected V E case, see Fig. (4). The 
chaotic motion corresponds to the trajectories that reach the non convex part 
of the boundary, including its internal part (Fig. (5)). The trajectories that 
are restricted to convex parts of the boundary are regular. 

The dependency of our results on a shows no surprises. The smaller is the 
parameter a, the close to E c we have the onset of chaos. For a given a, the 
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Fig. 3. Some trajectories in the configuration space corresponding to the Poincare's 
section of Fig. (2): (a), (b) and (c) the regular ones, corresponding to trajecto- 
ries that do not probe the non convex part of the boundary; (d) a chaotic one, 
corresponding to a trajectory that bounce on the non convex part of the boundary. 




-0.5 0.5 1 1.5 2 2.5 3 3.5 



Fig. 4. Typical Poincare's section (x,p x ) across the plane y = for the multi- 
ply-connected T>e case. For this section, one has a = 1/2 and E/Eq = 1.02. 

value of E / Eq marking the onset of chaos can be predicted analytically. Since 
the chaotic motions arise from the bouncing on the non convex part of the 
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Fig. 5. Some trajectories in the configuration space corresponding to the Poincare's 
section of Fig. (4): (a) and (b) the regular ones, corresponding to trajectories that 
do not probe the non convex parts of the boundary; (c) and (d) the chaotic ones, 
trajectories bouncing on the non convex parts. 

boundary, the onset of chaos will be determined by the appearance of the non 
convex region, which, as one can easily see, is related to the intersections of 
the curves V(x,y) = E and x 2 + y 2 = a. These curves have two intersections 
in the non convex case, and no intersections in the convex one. The value of 
E corresponding to only one intersection marks the transition. 

The situation for the case of two or more monopoles is similar to the expounded 
here. Chaotic behavior always appears for the trajectories that probe the non 
convex part of &D E , including the internal parts of the boundary for the 
multiply-connected case. 



4 Conclusion 



We have shown that the local stability criterion based on the geodesic devi- 
ation equation for the Jacobi metric is neither necessary nor sufficient to the 
occurrence of chaos. For this purpose, we introduce a class of 2-dimensional 
systems exhibiting chaotic behavior, confirmed by numerical evaluation of 
Poincare sections, with everywhere positive Gaussian curvature, and thus lo- 
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cally stable. The chaotic behavior arises for certain trajectories bouncing on 
non convex parts of the boundary of the effective Riemannian manifold, re- 
sembling the behavior of some classical billiards [17]. Despite the fact that 
the positivity of the Gaussian curvature K ensures the convergence of nearby 
geodesies, chaotic behavior can arise due to defocusing of geodesies bouncing 
on non convex parts of the boundary 8T>e- Our result questions, once more, 
the viability of local, curvature-based criteria to predict chaotic behavior. 
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